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Abstract 

We present a fast numerical technique for calculating the series impedance 
matrix of systems with round conductors. The method is based on a sur- 
face admittance operator in combination with the method of moments 
and it accurately predicts both skin and proximity effects. Application 
to a three-phase armored cable with wire screens demonstrates a speed- 
up by a factor of about 100 compared to a finite elements computation. 
The inclusion of proximity effect in combination with the high efficiency 
makes the new method very attractive for cable modeling within EMTP- 
type simulation tools. Currently, these tools can only take skin effect into 
account. 



1 Introduction and Motivation 

Electromagnetic transients have a significant impact on the design, operation 
and performance of electrical power systems [U [2] . Transients can be caused by 
several phenomena such as lightning discharges, breaker operations, faults, and 
the use of power electronics converters. Since electromagnetic transients involve 
a wide band of frequencies, ranging from DC up to the low MHz range, their 
simulation requires broadband and accurate models for all network components 
to predict the network response outside the nominal sinusoidal regime. For 
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the modeling of underground cables, it is necessary to calculate the per-unit- 
length (p.u.l.) cable series impedance matrix [3] over a wide band of discrete 
frequencies while taking into account frequency-dependent phenomena such as 
skin and proximity effects. The series impedance is next used as input data for 
alternative frequency-dependent cable models |1],[5]. 

Traditionally, the series impedance is calculated using analytic formulas 
which account for skin effect only, since they assume a symmetrical distribution 
of current in all conductors [3]. These formulas are combined with systematic 
procedures for computing the series impedance matrix for multi-conductor sys- 
tems [5]. Although simple and highly efficient, this traditional approach neglects 
proximity effects. Proximity-aware formulas are available only for two conduc- 
tor systems [3 13] . Ignoring proximity effect is acceptable for overhead lines and 
for widely spaced single core cable systems. However, in the case of three-phase 
cables, pipe-type cables and closely packed single core cables, the small distance 
between conductors in combination with the non-coaxial arrangement leads to 
significant proximity effects. This issue is relevant also in the modeling of um- 
bilical cables for offshore oil and gas power supply and control Ignoring 
proximity effect leads to an underestimate of the cable losses at the operating 
frequency, and the transient waveforms are also affected, in particular in sit- 
uations where waves propagate between the screens and between the screens 
and ground [51 HO] ■ The latter situation is highly relevant in the simulation of 
cross-bonded cable systems. 

Because of the limitations of analytic formulas, several numerical techniques 
have been proposed. The harmonic expansion method accounts for proxim- 
ity effects at high frequencies under the assumption that the skin effect is 
fully developed [31 [TTl [TH [T3] . Techniques based on the finite element method 
(FEM) [HI m] fully predict proximity effect at both low and high frequencies, 
but they tend to be excessively time-consuming because of the fine mesh re- 
quired to properly discretize the cross session. A similar issue arises with tech- 
niques based on conductor partitioning [15l [TBI [T71 [181 [T9| . Moreover, with both 
FEM and conductor partitioning techniques, the discretization must be refined 
as frequency increases to properly capture the pronounced skin effect, further 
reducing computational efficiency 

In this paper, we overcome these issues by developing a fast numerical tech- 
nique for computing the series impedance matrix of cables with round conduc- 
tors. Our approach, denoted henceforth as MoM-SO, combines the method of 
moments (MoM) with a surface admittance operator (SO) introduced in [2U] 
for rectangular conductors. The proposed approach is not available elsewhere 
in the literature, since the authors of [5^, after presenting the surface admit- 
tance operator for both rectangular and round conductors, focus their attention 
on the first case. Through the surface operator, we replace each conductor with 
the surrounding medium, while introducing an equivalent current density on 
the surface of each conductor. The current density and the longitudinal elec- 
tric field on all conductors are then related through the electric field integral 
equation |21) , and their spatial dependence is expressed in terms of Fourier com- 
ponents. Finally, the method of moments [22j is applied to compute the p.u.l. 
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parameters of the line. Numerical results will show that the proposed approach 
is much more efficient than state-of-the-art FEM techniques, since few Fourier 
components are sufficient to accurately model skin and proximity effect at any 
frequency. Moreover, the discretization of the electric field integral equation is 
entirely perfomed using analytic formulas, avoiding numerical integration used 
in previous works [131 HO] ■ This achievement further improves the robustness 
and speed of MoM-SO. 

The paper is organized as follows. In Sec. [5] we state the problem from 
a theoretical standpoint and review the two fundamental relations exploited in 
this work, namely the electric field integral equation and the surface admittance 
operator. In Sec. |3j the two relations are discretized with the method of mo- 
ments. In Sec. |4j the new MoM-SO method is first validated against analytical 
formulas for a simple two-conductor system, and then compared against FEM 
on a three-phase armored cable. 

2 Impedance Computation via a Surface Admit- 
tance Operator 

2.1 Problem Statement 

We consider a transmission line made by P round conductors parallel to the z 
axis and surrounded by a homogeneous medium. An example of line cross sec- 
tion is shown in Fig. [l] where the permittivity, permeability and conductivity of 
the conductors are denoted as e, /x, and a respectively. The surrounding medium 
is assumed to be lossless, with permittivity Sout and permeability /io- Our goal 
is to compute the per-unit-length (p.u.l.) resistance ^{lo) and inductance ^{lo) 
matrices of the line as defined by the Telegraphers' equation [3] 

1^ ^-[^H + j...5fM]i (1) 

az 

where V — \Vi - ■ ■ Vp]^ is a P x 1 vector collecting the potential Vp of each 
conductor. Similarly, the P x 1 vector I = [/i • • • Ip]^ is formed by the current 
Ip in each conductor. The parameters S^{lS) and JSf (w) in ([T]) are commonly 
referred as partial p.u.l. parameters [23] . From them, one can easily obtain 
the p.u.l. parameters of the line with any conductor taken as reference for the 
voltages and as return path for the currents 

2.2 Surface Admittance Operator 

In order to compute the p.u.l. impedance of the line, we follow the approach 
of [20] which relies on a surface admittance operator. We replace each conductor 
with the surrounding medium and, to maintain the electric field Ez outside the 
conductors' volume unchanged, we introduce a surface current density Js on 
their contour, as shown in Fig. [2| The electric field inside the conductors' 
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Figure 1: Sample line cross section. The conductors are depicted in white 
(permittivity: e, permeability: /i, conductivity: a). The surrounding medium 
is depicted in grey (permittivity: £out^ permeability: /ip)- The unit vector fi is 
the normal to the conductors surface. 




Figure 2: Line cross section after application of the equivalence theorem. The 
conductors medium is replaced by the surrounding medium, and an equivalent 
current density is introduced on the surface of the conductors. The contour 
and radius of the p-th conductor are denoted with Cp and Cp, respectively. 



volume takes instead a fictitious value Ez. The current density Jg, directed 
along z, can be found with the equivalence theorem ^ and reads 



J a ~ Hf. — Ht 



(2) 



where Hf is the component of the magnetic field tangential to the conductor's 
surface, evaluated before the application of the equivalence theorem (configura- 
tion of Fig.jlJ. Ht is the same quantity evaluated after the equivalence theorem 
has been applied (configuration of Fig. [2| . 

On the boundary of each conductor, the tangential component of the mag- 
netic field is related to the longitudinal electric field [20] 



1 dE, 
jujfi dn 



(3) 



where ^ denotes the directional derivative 



with respect to the unit vector 
n normal to the conductors surface. Similarly, for Ht we can write 



1 dE, 



juj^o dn 

By substituting ^ and Q into we obtain the relation 



fi dn 



1 dE, 
^0 9n 



(4) 



(5) 
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which defines a surface admittance operator that maps the electric field on the 
conductor's boundary onto the equivalent current density Js- Equation (|5| 
extends the formula given in [20] to the case of magnetic conductors ^ /xq). 



2.3 Electric Field Integral Equation 

After applying the equivalence theorem to each conductor, the medium becomes 
homogeneous, and we can easily relate the equivalent current density Jg to the 
electric field by means of the electric field integral equation [21] . This process 
involves only the boundary Cp of each conductor, that we describe with the 
position vector 

rp{6) — {xp + ap cos 9)x + {yp + ap sin 9) y , (6) 

where is the azimuthal coordinate, (xp, yp) are the coordinates of the conductor 
center, and Op is the conductor radius, as illustrated in Fig. [2[ The unit vectors 
X and y are aligned with the x and y axis, respectively. We denote the equivalent 
current density and the electric field on the boundary of the p-th conductor as 
Js^\9) and E'f\o), respectively. Using the electric field integral equation, we 
can express the electric field E^^^ (6) on the surface of the p-th conductor as 

EiP\0) = E / Ji'\0')G {fp(e),r,{0')) a,d9' 

dV 
oz 

where the first term accounts for the field generated by the current on each 
conductor, while the second term is related to the scalar potential V. The 
integral kernel 

G(fp, fq) In \fp - fq\ (8) 

is the Green's function of an infinite space [53] . On the p-th conductor, the scalar 
potential V is equal to the conductor potential Vp that appears in the Teleg- 
raphers' equation ([T]). Therefore, we can replace the last term in ^ with ([!]), 
obtaining 



^ r 



^\ji''\e')G{rp{e),T,{e'))a,de' 



p 

+ [^p,(c^) + j^^pqiu)] Iq , (9) 

9=1 

where 3^pq{ui) and Jfpq{uj) are the elements in position {p,q) of the matrices 
Mito) and jSf (w), respectively. Equation ^ combined with the surface admit- 
tance operator ([5| will allow us to compute the p.u.l. parameters of the line. 
In the next Section, we introduce a discretization of these relations suitable for 
numerical computations. 
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3 Numerical Formulation 



3.1 Discretization of the Surface Admittance Operator 

Given the cylindrical geometry of the conductors, we approximate the field and 
the current on the p-th conductor by means of a truncated Fourier series 



E 



(p) 



(9) - E 



(p) gjne 



1 

27ra„ 



E 



(10) 



(11) 



The truncation order Np controls the accuracy and the computational cost of 
the numerical technique. Examples will show that a low Np, of the order of 
2 — 3, delivers very accurate results while minimizing the computation time. 
Owing to the normalization factor , the total current Ip flowing in the p-th. 
conductor is simply given by the constant term of the series |20j 



ip) 



(12) 



(v) (v) 

When E'/' and Jr' are expressed in Fourier series, the surface admittance 



operator ([S]) can be rewritten in terms of the Fourier coefficients as 



j(p) = ^(P) 



27r 



(13) 



of order \n\, and J|'„|( 



where J|„| (.) is the Bessel function of the first kind 

is its derivative. The wavenumber in the conductors and in the surrounding 
medium are given respectively by 



k = ■\/ uj^,[uje — ja) 

out • 



(14) 
(15) 



In order to simplify the oncoming equations, we introduce a compact matrix 



(v) 

notation. We collect all field coefhcients En in the column vector 

T 



E 



E 



(1) 



E 



(1) 

Ni 



E 



(2) 

-N2 



E 



(2) 

N2 



(16) 



(v) 

and all current coefficients J„ in 



J = 



^(1) 



j(2) 



The vectors E and J have size 



-N2 



N = 



p 

E 



J 



(2) 



(2^p + 1) , 



(17) 



(18) 
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the total number of field and current coeffidents. Relation (13) can be written 
in terms of E and J as 

J = Y,,E . (19) 

where is a diagonal matrix. This matrix is the discrete version of the surface 
admittance operator defined by Q for each round conductor. Finally, ( 12 ) can 
be written in terms of I and J as 

I = U"^ J , (20) 

where U is a constant N x P matrix made by all zeros and a single "1" in each 
column. In column p, the "1" is in the same row as the coefficient Jq^-* in (17). 

3.2 Discretization of the Electric Field Integral Equation 

We now cast the electric field integral equation (|9| into a set of algebraic equa- 
tions using the method of moments, a numerical method to solve integral and 
differential equations [22]. First, we substitute (10) and (11) into (|9|, obtaining 



E Ji'' r e=-''G{r,{e),r,{B'))d9' 



e 



P 

E [^pqH + J^-Sfp,(u;)] Iq . (21) 

q=l 



Then, we project (21) onto the Fourier basis functions e-*" by applying the 
operator 

r27r 

[.]e-^''''^de n' = -Np, ...,Np (22) 
to both sides of the equation, obtaining 

^!?^ = ^-MoE E G2;;f4')+ 

P 

Y,[^p,{u)+juj^p^{Lo)]Iq, (23) 



for p = 1, . . . , P, where 



I 1 when n' = 
I when n , 
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and where G^^^'' denotes the {n',n) entry of the matrix G(P'«l This matrix 
describes the contribution of the current on the q-th conductor to the field on 
the p-th conductor. The entries of G'P''^^ are given by the double integral 



2-rr /.27r 



G (fp(0), f,(0')) eJ'("^'-"'^)d^d0' , (25) 



which can be computed analytically as shown in the Appendix. Using the matrix 
notation set in ( 16 1 and ( 17 1 , the equations ( 23 1 can be written in compact form 
as 

E = jcj/ioGJ + U + jiu^iuj)] I , (26) 



where G is the block matrix 



G 



G(^'i) 



(27) 



Equation (26) is the discrete counterpart of the electric field integral equa- 
tion (|9|. 

3.3 Computation of the Per-Unit-Length Parameters 

The p.u.l. parameters of the line can be obtained by combining the discretized 
surface admittance operator ( 19 ) with the discretized electric field integral equa- 
tion ( 26 ) as follows. First, we left multiply ( 26 ) by Yg and, using ( 19 1, we obtain 

(28) 



J = jWyUoY^GJ + Y,U [M{uj) + jwJf (cj)] I . 
The current density coefficients J can be expressed as 

J = (1 - jcj^oYsG)-iY,U [^{uj) + jwJSf (cj)] I , 



(29) 



where 1 denotes the N x N identity matrix. Left multiplication of ( 29 1 by U"^ 
leads to the following expression for the line currents I 



(30) 



Since (30) must hold for any I, the product of the two expressions inside the 
square brackets must be equal to the identity matrix. Consequently, we have 
that 

^(w) + jwjSf (w) = [UT(1 - juj^nYsG)-^Ys\5] . (31) 



By taking the real and imaginary part of (31 ) we finally obtain the formulas for 
computing the p.u.l. resistance and inductance matrices 

^(^) =Me{[uT(l-jc^^oY,G)-iY,U]"'} , 



UJ] 



-ilm{[uT(l-j^^oY,G)-^Y,U] '} 



(32) 
(33) 
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4 Numerical results 



4.1 Two Round Conductors 

In order to validate the proposed technique against analytic formulas, we con- 
sider a line made by two parallel round conductors with radius a = 10 mm made 
of copper (cr = 58 • 10^ S/m, ^ = /xg). Two different values for the center-to- 
center distance between the wires have been used, namely D = 100 mm and 
D = 25 mm. In the first case proximity effect is negligible, due to the wide 
separation. In the second case it is instead significant. 

The p.u.l. computed with a MATLAB implementation of MoM-SO have 
been compared against two different sets of analytic formulas. The first set is 
valid at high frequency [3] because it assumes a fully-developed skin effect, and 
gives the p.u.l. resistance and inductance as 

(34) 

^,,, = /^cosh-i(^£) , (35) 
where Rg — (crS)^^ is the surface resistance and 

^ (36) 



is the skin depth. The second set of formulas accounts for the frequency de- 
pendence of the internal impedance of the wire 3fint , which can be calculated 
analytically under the assumption of wide separation [3] 

= ^ ber(e)+jbei(0 ^ ^^^^ 
V^naaS bei'($) - j ber'(^) ' 

where ^ = V^f- The Kelvin functions ber(^) and bei(^) are the real and 
imaginary part of J^o{£,e^^^), respectively [35]. The total p.u.l. impedance of 
the line is thus 

^ = 22fint + j<-^-^ext ■ (38) 



in (371 



Formulas (34)-(35) account for proximity effect, which is instead neglected 



The p.u.l. parameters have been computed from 1 Hz to 1 MHz with the 
truncation order Np set to 3. No noticeable changes have been observed beyond 
this value. The computation of the parameters took 15 ms per frequency sample 
on a 3.4 GHz CPU. The discretization of the Green's function which leads 
to the matrix G took less than 5 ms. Figures [3] and |4] compare the p.u.l. 
parameters computed with MoM-SO against the results obtained from the two 
analytic formulas. In all cases, the numerical results correctly approach the 



exact high frequency value given by formulas ( 34 ) and ( 35 ) . In the case of wide 
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Table 1: Characteristics of the cables in the example of Sec. 4.2 



Item 


Parameters 


Core 


(T = 58 • 10^ S/m, r = 


10.0 


mm 


Insulation 


t — 4.0 mm, Ej, = 


2.3 




Wire screen 


32 wires, r = 0.5 mm, a — 


58 • 


10^ S/m 


Jacket 


t = 2 mm, er = 


2.3 





Table 2: Arm or characteristics for the structure conside red in Sec. 14.21 



Item 


Parameters 


Armor outer diameter 


88.26 mm 


Wire diameter 


3 mm 


Conductivity 


10^ S/m 


fir 


100 


N.o. wires per layer 


70 



separation, shown in the top panels, the numerical results also correctly predict 



the frequency-dependent behavior of the wire's internal impedance (37) due to 
skin effect. A little discrepancy is visible in the inductance at low frequency 
(top panel of Fig. [4]), due to the small but not negligible proximity effect. The 
error introduced by p7[ ) becomes more significant when the wires separation 
is reduced to 25 mm, as shown in the bottom panel of Figures [3] and |4j The 
non-uniform current distribution induced by the wires proximity is visible in 
Fig. [sj which also shows the development of skin effect. 

4.2 Three-Phase Armored Cable 

We consider the three-phase armored cable in Fig. [6] which features three wire 
screens and a steel armoring, for a total of 293 circular subconductors. The 
key parameters are listed in Tables [T] and [2| respectively. Using MoM-SO we 
computed the 3x3 series impedance matrix with respect to the three phase 
conductors with the screens continuously grounded along the cable. 

Table |3] shows the calculated positive and zero sequence resistance and re- 
actance per km. The computation has been performed with three different 
truncation orders: Np = 0, Np = 3 and Np = 7. As a validation we used a FEM 
implementation [81 with a very fine mesh (177,456 triangles). It is observed that 
with orders Np — 3 and Np — 7 we get a result which deviates by less than 1% 
from the FEM result. 

Figs. [7] and [8] show the positive and zero sequence resistance and inductance 
as a function of frequency, from 1 Hz to 1 MHz. It is observed that with Np = 0, 
significant errors result as the proximity effects are ignored. Indeed, by setting 



Np = in ( 10 ) and ( 11 ), one assumes a circularly-symmetric current distribution 



10 



Table 3: Positive- and zero-sequence impedance of the three-phase cable of 
Sec. 4.2 at 50 Hz. MoM-SO is compared against a finite element approach [8]. 





MoM-SO (proposed) 


FEM 


Nj, = 


Np = 'S 


Np = 7 


i?+[0/km] 


0.06905 


0.07259 


0.07261 


0.07218 


X+[n/km] 


0.08703 


0.09042 


0.09048 


0.09041 


Ro[n/km] 


0.2386 


0.2437 


0.2438 


0.2459 


Xo[fi/km] 


0.08033 


0.08943 


0.08958 


0.08975 



on the conductors. With Np = 3 and Np — 7, a. virtually identical result is 
achieved which agrees very well with the FEM result. At very high frequencies, 
however, the FEM result deviates somewhat from that of the proposed approach 
since the mesh division is not sufficiently fine to properly account for the very 
small skin depth. In the proposed technique, instead, skin effect is implicitly 
and fully described by the surface admittance operator, and does not affect 
the discretization of the problem, which depends only on the proximity of the 
conductors. As a result, the level of discretization, controlled by Np, does not 
have to be increased as frequency grows, making MoM-SO much more efficient 
than FEM. 

Timing results, reported in Table |4j demonstrate the excellent performance 
of the developed algorithm. With MoM-SO, there is first a computation time 
for the Green's function discretization of 11.6 s {Np = 3) or 16.5 s [Np — 7). 
Then, for computing each frequency sample one needs 2.01 s for Np = 3 or 15.5 s 
for Np = 7. Since Np — 3 was found sufficient for obtaining accurate results, 
the total computational cost for computing the 31 samples in this example is 
T = 11.6 31 X 2.01 = 73.9 s. The computation time using FEM is much 
higher, requiring 440 s per frequency sample. This is 220 times slower than the 
per-sample computation time of 2.01 s using the new method with Np = 3. We 
can therefore safely state that MoM-SO is at least 100 times faster than the 
FEM approach when several frequency samples are needed. 



5 Discussion 

The new method is applicable for any system of non-overlapping solid, round 
conductors. In practice, cables often include tubular conductors (sheaths) as 
well as ground return in earth. The authors are currently extending the method 
to handle also these situations. As for the example of the three-phase cable, the 
modeling in this paper is fully applicable since the sheath conductor consists of 
wires and because the thick armor permits to ignore the external medium. 
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Table 4: Three-phase cable example of Sec. 4.2 computation time for MoM-SO 
and FEM. 





MoM-SO (proposed) 


FEM 


Np = 


Np = 3 


Np = 7 


Green's 
function dis- 
cretization 


11.6 s 


13.7 s 


16.5 s 




Impedance 
computa- 
tion (per 
frequency 
sample) 


0.085 s 


2.01 s 


15.5 s 


440 s* 



All computations were performed on a system 
with a 2.5 GHz CPU and 16 GB of memory. 
*Mesh size: 177,456 triangles. 



6 Conclusion 

We presented a very fast and accurate algorithm for computing the series 
impedance of systems of round conductors. The method combines a surface 
operator with the method of moments which permits to compute the complete 
series impedance matrix while taking into account both skin and proximity ef- 
fects. This capability is of major importance in cable system modeling as the 
short lateral distance between cables often leads to significant proximity effects. 
Due to its efficient discretization of the underlying electromagnetic problem, the 
algorithm outperforms state-of-the-art techniques based on finite elements by a 
factor of about 100. The computed resistance and inductance can be used for 
an accurate prediction of electromagnetic transients in EMTP-type programs 
when combined with an appropriate frequency-dependent cable model. 



Analytical evaluation of the Green's matrix G 

The discretization of the Green's function requires the computation of the double 



integral ( 25 ) . After substitution of ([8|) , the integral reads 

r27r 



1 



(2^)2 



fn{e)e 



'jn'e 



de. 







where 



fn{e) 



2tt 



27r 



(39) 
(40) 
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We first calculate the integral in ( 40 ) , which can be expanded using ^ to obtain 



fn{0) 



2-K 



In 



1 

An 



2ir 



In [(r")' + al- 2a,r" cos(0" - 9')] e''''^' d9' , (41) 



where r"{6) is the auxiliary vector 

r"(6l) = rp{e) - [xqx + yqy) = r"{cose"x + &me"y) , 



(42) 



which is constant with respect to the integration variable 0' . We denote its 
modulus and its angle with r" and 9" respectiveljQ The solution to the last 
integral in ( 41 ) is presented in [3] and reads 



fn{e) 



ln(r") 



2\n\(r" 



for n — Q , 
for n ^ . 



(43) 



Next, we solve the integral in (39). The solution involves several cases, which 
are itemize for better readability: 



p ^ q , n = 0: in this case, the integral in (39) is analogous to (40) and 
can be solved with the formulas given in ,3,. When n' = 0, we have 



1 



(44) 



and when n' 7^ we have 



1 



Att In' 



where Xp,, = {xp - Xq), yp^q = [yp - yq), and 



(45) 



(46) 



p q , n > 0, n' > 1: if we manipulate the integrand function as follows 



fni9) 



(a,)"eJ"^ 
2n(r")" 

2n 



2n 



r e 



(r")2 



Xp^q + ttp COS 6 + j{yp,q + ap sin t 



{xp^q + ttp COS + [yp^q + Op sin Oy 



1 



2n {Xp^q - jyp^q + ttpC J")" ' 



^For the sake of clarity of the notation, we omit from r" and d" the dependence on 1 
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we can rewrite ( 39 ) as a complex integral 



G 



dz . 



(47) 



performed over the closed path z — apt with 9 = [0,27r]. Since we 
assume n' > 1, the integrand function in (47) will have no poles inside the 



integration path and, using the residue theorem 27J, we have that 



(48) 



p 7^ g , n > 0, n' < 1: when n' < 1, the integrand function in (47) has a 



pole in z = —Xp^q + jyp^q- Applying again the residue theorem we obtain 

(49) 



^(p,,) ^ -niaqr (n-n'- l\ [d, - jd,)-"+" 



with (") being the binomial coefficient. 



(27r)2n 



• p ^ q , n < 0: the Green's matrix entries for n < can be obtained 
from ( 48 ) and ( 49 ) by symmetry 



(50) 



where * denotes complex conjugation. This identity follows from the sym- 
metry relation /-„(6l) = /„(6')*. 



• p — q: in this case we have that Xp^q = yp^q — 0, and using (47) one can 
easily show that 

i(p,p) _ /^Inap ifn' = 0, 
if 7i' ^ , 



q(p,p) _ 
n' ,n 



4^ I n I 



if 71 7^ 0, n' = n , 
if n 7^ 0, n' 7^ n . 



(51) 
(52) 
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Figure 3: P.u.l. resistance of the two wires line of Sec. |4.1| computed with 



MoM-SO (crosses) , the analytic formula ( 34 ) valid at high frequency (dash-dot 



line), and formula (38) (solid line). Wires separation is 100 mm (top panel) and 



25 mm (bottom panel). 
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Figure 4: P.u.l. inductance of the two wires line of Sec. |4.1| computed with 



MoM-SO (crosses) , the analytic formula ( 35 ) valid at high frequency (dash-dot 



line), and formula (38) (solid line). Wires separation is 100 mm (top panel) and 



25 mm (bottom panel). 
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(a) 515 Hz 



(b) 3.63 kHz 



Figure 5: Two wires example of Sec. 4.1 in the case of a separation of 25 mm: 
current density plot for two different frequencies, obtained with MoM-SO. 
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Figure 7: P.u.l. resistance of the three-phase cable of Sec. |4.2[ obtained with 
MoM-SO and FEM. For MoM-SO, three different truncation orders Np are 
considered. 
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Figure 8: P.u.l. inductance of the three-phase cable of Sec. |4.2[ obtained with 
MoM-SO and FEM. For MoM-SO, three different truncation orders Np are 
considered. 
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